Accurate determination of inflationary perturbations 
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We use a numerical code for accurate computation of the amplitude of linear density perturbations 
and gravitational waves generated by single-field inflation models to study the accuracy of existing 
analytic results based on the slow-roll approximation. We use our code to calculate the coefficient 
of an expansion about the exact analytic result for power-law inflation; this generates a fitting 
function which can be applied to all inflationary models to obtain extremely accurate results. In 
the appropriate limit our results confirm the Stewart-Lyth analytic second-order calculation, and 
we find that their results are very accurate for inflationary models favoured by current observational 
constraints. 
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I. INTRODUCTION 

In anticipation of a near-future launch of a satellite 
capable of measuring microwave background anisotropics 
to an accuracy of a few percent or better across a wide 
range of angular scales [IJ, attention has recently been 
directed towards obtaining highly accurate predictions 
of the anisotropies for given cosmological models. It is 
now possible for fast and extremely accurate calculations, 
good to a percent or so, of the radiation anisotropy power 
spectrum to be made , given the spectrum of density 
perturbations in the universe. 

A very promising theory for the generation of such 
perturbations is cosmological inflation ||, wherein the 
universe experiences a period of accelerated expansion 
in its very early stages. In order to take advantage of 
the great observational accuracy anticipated, it is desir- 
able to make as accurate calculations as possible of the 
perturbations which inflation produces. Such accurate 
calculations are the focus of this paper. 

When density perturbations from inflation were first 
discussed , the result was typically quoted along the 
lines of 
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where the right hand side is evaluated as the scale crosses 
outside the Hubble radius during inflation, and the left 
hand side as it crosses back inside after inflation. Here H 
is the Hubble parameter and </> is the scalar field driving 
inflation. That the result was given only approximately 
was through a combination of the calculations not being 
completely accurate, and the lack of a precise definition 
of the left hand side of this equation, usually described 
as the fractional density perturbation when a scale enters 
the horizon. Fairly quickly, however, a precise definition 
of this quantity was made and the coefficient computed 
^,[1. In more modern notation, defined later in this pa- 
per, the result is 



The right hand side is evaluated as before, but now the 
left hand side, giving the spectrum V of the curvature 
perturbation 1Z, holds at any time when the scale is much 
larger than the Hubble radius, the curvature perturba- 
tion 1Z being constant in that limit. The uncertainty in 
the calculation was then due to the fact that its deriva- 
tion depended on the inflationary slow-roll approxima- 
tion; for some models this can be extraordinarily good 
but in many the associated error would be expected to 
be order of ten percent or greater. 

The next development was the presentation of a pre- 
cise calculation in the case of power-law inflation, by Lyth 
and Stewart |?J . The analogous but simpler exact calcu- 
lation for the gravitational wave amplitude in this model 
was made much earlier by Abbott and Wise ||. The 
precise calculation allowed a direct evaluation of the un- 
certainty induced in these models by use of the slow-roll 
approximation. Stewart and Lyth then went on to use 
this exact result to analytically compute the next-order 
slow- roll correction to the standard formula || . In princi- 
ple, one would like to carry out the expansion about any 
power-law inflation model, but in practice analytic con- 
siderations required a restriction to an expansion about 
power-law inflation models which were themselves near 
the de Sitter limit. One should therefore regard that 
calculation as a slow-roll expansion rather than an ex- 
pansion about the power-law inflation solution. 

A slightly different kind of calculation, which we shall 
not further discuss, can be done in the case where the po- 
tential takes the form V(<j>) = Vq ± m 2 (j) 2 /2, in the limit 
where Vq dominates. For the inverted harmonic oscilla- 
tor (natural inflation) this was carried out by Stewart and 
Lyth Q and for positive mass (hybrid inflation) an anal- 
ysis was performed by Garcfa-Bellido and Wands Jl0| . 
There a much more accurate calculation can be made 
than the slow-roll one, even if the mass is quite large, 
but again this calculation is restricted to this particular 
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class of models. 

For most models of inflation, the Stewart -Lyth cal- 
culation should give extremely accurate predictions. In 
this paper, we use a numerical code to test their accu- 
racy against exact results, finding that the agreement is 
extremely good. We then use our code as the basis for 
an even more accurate scheme, using it to carry out an 
expansion about any power-law inflation solution, even 
those far from the de Sitter limit. This gives a highly 
accurate fitting function for the amplitude in any infla- 
tion model, by regarding them as expansions about the 
'nearest' power-law inflation model. 



II. THE DENSITY PERTURBATION 
CALCULATION 

A. Evolution of perturbations 

Standard calculations of the density perturbation am- 
plitude from inflationary models rely on two separate 
types of assumption, linear perturbation theory and the 
slow-roll approximation for the inflationary evolution of 
the background space-time. The first of these approxi- 
mations is considerably better than the second, since we 
know from COBE that the amplitude of perturbations is 
about one part in 10 5 , and we shall use it throughout. 

The slow-roll approximation is typically much less 
good, because in order for inflation to end the slow-roll 
conditions must break. While for some inflationary mod- 
els it is extremely good, it is perfectly possible for it to fail 
by ten percent or more in others, including some still per- 
mitted by present observational constraints. Hence our 
aim in this paper to make calculations which do not de- 
pend on the slow-roll approximation, and also to provide 
more accurate approximations for the density perturba- 
tions based on a version of the slow-roll expansion. 

Before proceeding to the calculation, let us highlight 
that we are limiting discussion to inflation driven by a 
single, canonically normalized, scalar field in Einstein 
gravity. While this looks like quite a restrictive assump- 
tion, is is actually fairly general, in the sense that almost 
any extended theory of gravity, such as a higher-order or 
scalar-tensor theory, can be rewritten as Einstein grav- 
ity via a conformal transformation. Further, models with 
more than one scalar field often effectively only have a 
single degree of freedom, with fluctuations transverse to 
the evolving field having their fluctuations suppressed by 
a large effective mass; this situation occurs, for exam- 
ple, when extended inflation is conformally transformed 
to give power-law inflation jO], and also in the hybrid 
inflation model [|l2fl . However, our calculations will not 
be valid in any theory in which fluctuations in more than 
one field are important, and it is possible to construct 
situations where this is the case [|13[ . 

Our analysis is based on the formalism devised by 
Stewart and Lyth H] to carry out a second-order analytic 



calculation. We shall not reproduce the derivation of the 
equations here, instead simply reproducing the important 
ones.0 With the usual notation of a for scale factor, H 
for Hubble parameter, ip for the inflationary scalar field 
and overdot as derivative with respect to cosmic time, 
they introduce a new quantity z defined by 
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The quantity one desires to calculate is the curvature 
perturbation 1Z, defined as in Ref. ||. It is convenient to 
define a related quantity u, defined by 



u = —zTZ , 
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which is a gauge invariant potential |p_ 5|| . 

A Fourier expansion of u into comoving modes itk is 
carried out, and these can be_ shown to obey the remark- 
ably simple equation JlqJlF 
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where r is the conformal time defined by dr = dt/a{t) 
and k = |k| is the modulus of the wavenumber. Defining 
the spectrum V-nik) of the curvature perturbation in the 
standard way as 
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which statistical isotropy permits to depend only on the 
modulus k. From here on, we'll label modes by their 
modulus k rather than their full wavenumber. 

During inflation, a comoving scale evolves from well 
inside the Hubble radius to well outside it. Canonical 
quantization of u demands that in the initial state, where 
the expansion can be ignored, it approaches the standard 
flat space result H 
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In the opposite regime, where k can be ignored in Eq. (|^), 
the equation can be integrated directly to give 



Uk oc z , 
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implying that the spectrum Viz(k) approaches a constant 
in this regime. 



*An expanded discussion of the Stewart-Lyth calculation is 
given in Ref. juj. 
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Notice that the entire quantum part of this calculation 
is in the normalization at large k, which is the Minkowski 
space limit. Once that has been determined, then in a 
Heisenberg picture the state vector is time independent 
and the mode functions Uk obey the classical equations 
of motion which allow us to evolve them. The full theory 
of quantum fields in curved space-time is not required. 

The crucial question, of course, is the value of the pro- 
portionality constant in Eq. (||). We have seen that in 
both the short and long wavelength regimes the evolution 
of amplitude is independent of the way in which the uni- 
verse expands. The crucial feature which determines the 
final amplitude is therefore the way in which the universe 
evolves as the mode switches from one asymptotic regime 
to the other. This requires a modelling of the relevant 
inflationary epoch. 



B. Inflation and the slow-roll approximation 

Given a specific inflationary potential V{<j>) for the 
scalar field <fi, there is a well posed problem for the evo- 
lution of the background space time. The relevant be- 
haviour that we want to extract is that of (d 2 z/dr 2 )/z, 
as required in Eq. (||). 

The most compact way of describing inflationary solu- 
tions is the Hamilton-Jacobi formalism flTj] , where one 
adopts the scalar field itself as the time variable and in 
which the solution is described via H((f>). The equations 
of motion read 
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If required, the potential corresponding to any solution 
H(4>) is easily obtained from the first of these. 

From this fundamental quantity, we can introduce a se- 
ries of functions containing higher and higher derivatives 
of H{4>) viaQ 



<<f>) 
77O) 



= '"pi 



EL 

IT 



H" 



4tt 

47T 



'"-PI 
47T 



H 



' H'H' 1 



1/2 



H 2 

if™ 



l/n 



(12) 

(13) 
(14) 

(n > 3) , (15) 



^The functions rj and £ are equivalent to /3i and fa; we give 
them a special name as they crop up frequently below. The 
first function, e, doesn't fit into the pattern of the rest. 



where prime is a derivative with respect to (j> and '(n+1)' 
symbolizes the taking of n + 1 derivatives with respect to 
(/>. These functions are known as the slow-roll parameters, 
introduced by Liddle, Parsons and Barrow Jjql . They 
form the basis for the slow-roll expansion, which allows 
arbitrarily accurate solutions to the dynamical equations 
governing inflation to be obtained. For our purposes, 
they are useful because we can replace the general func- 
tion H{4>) by the values of H and these slow-roll parame- 
ters at a single value of <p — this is equivalent information 
because, modulo questions of convergence, the full func- 
tion H((f) can be reconstructed via a Taylor expansion 
about this point. Because we are interested primarily in 
the behaviour as a scale k crosses the Hubble radius, we 
shall choose the value of <j> when k = aH. Notice that 
the calculation we carry out is for a single value of k; 
even with the same inflationary potential, the values of 
these expansion coefficients change if we change the scale 
k, since that corresponds to moving to a different part of 
the potential. 

As we mentioned at the start of the Section, the in- 
flationary input into the perturbation equation is that it 
determines (d 2 z / dr 2 ) / z in Eq. (0). This quantity can be 
rewritten in terms of the slow-roll parameters, as 
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Although this looks like it might be the start of an ex- 
pansion, it is actually exact. 



C. Exact solution for power- law inflation 

We begin by rederiving the exact solution describing 
power-law inflation, following Lyth and Stewart 
The evolution of the background space-time is given by 



H{(f) = H exp ( \ — 
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where Ho and p > 1 are constants. The scalar field is 
translated so that <j) = <p a at the time when the scale k we 
are interested in obeys k = aH. This case is the simplest 
because all the slow-roll parameters are constant, and in 
fact they are all equal to 1/p. This means that Eq. jl^ ) 
takes on a particularly simple form. 

We also require an expression for the conformal time, 
which for power-law inflation can be derived using a trick 
of integrating by parts 



dt 
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~a 2 H 
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which for constant e implies 
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With these results, Eq. (|5|) for the perturbations re- 
duces to a Bessel equation 



dr 2 



{v 2 ~ 1/4) 



u k = , 



where 
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is a constant. The solution with the correct short-scale 
behaviour, shown in Eq. (||), is 

u k (r) = ^ e W/2W2 ( _ r) i/ 2i/ (i) ( _ fcr) (22) 

where Hu is the Hankel function of the first kind of 
order v. 

The result we desire is the asymptotic form of the so- 
lution; taking k/aH — > gives the asymptotic form 
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where T is the usual gamma function. On substitution 
into the expression for the power spectrum, Eq. (Q), this 
gives 



' n (k) ~ Z r(3/2)^ m% x \H'\ 
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It should be stressed that, despite appearances, this equa- 
tion does not give the value of the perturbation as it 
crosses the Hubble radius. Instead, it gives the asymp- 
totic value as k/aH — ► 0, rewritten in terms of the values 
which quantities had at Hubble radius crossing. 

In the limit of small e, this approaches the standard 
result 
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Note however that the amplitude diverges in that limit. 



D. Strategy for a general calculation 

The aim is to use this exact solution to estimate the 
amplitude expected in any inflationary model, by relat- 
ing the expansion behaviour crucial for the perturbation 
generation on a given scale k to that of a suitable power- 
law inflation model. Power-law inflation in effect gives a 
two parameter set of solutions, the input parameters be- 
ing Hq and e — 1/p. In a general inflation model H and e 
vary in some way with time, and take on some particular 
values when k — aH (at <f) — cf>o). We can then choose a 
reference power-law inflation model which matches these 



values; this guarantees that the amplitude and slope of 
H(cf>) about (/)q are the same in the reference power-law 
model as in the true model. The perturbation amplitude 
is then computed in the reference model. Since the final 
answer is primarily determined by the behaviour when 
k ~ aH, this generates the appropriate approximate an- 
swer, provided the true H{<p) respects this linear approx- 
imation sufficiently accurately across the relevant scales. 
This reasoning leads to the standard result, Eq. (p5|). 

Notice that the calculation is carried out for a single 
scale fc; to generate the complete spectrum, then at each 
k there is a different reference power-law inflation model 
because in general both H and e evolve in a way different 
to power-law inflation. 

If one follows the logic of the preceding paragraphs, 
then it comes as a surprise that the standard result 
quoted is the small e result, Eq. j25|), rather than the 
exact power-law inflation result Eq. (24). One reason 



for this is that the level of accuracy of including the e- 
dependent prefactor has usually not been required. How- 
ever, a more important reason is that the prefactor in 
Eq. (^4|), which depends only fairly weakly on e, is not 
general enough; there are terms connected to r\ (i.e. from 
the second derivative of H) which give comparable cor- 
rections. For a more accurate result than Eq. (63), these 
need to be taken into account. 

Stewart and Lyth ^| used the power-law inflation so- 
lution as the basis for an analytic expansion intended to 
hold for all inflation models; in an approximation where 
e and r\ are treated as negligibly varying they obtained 
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where C = — 2 + In 2 + 7 ~ —0.73 is a numerical constant 
(7 being the Euler constant). The slow- roll parameters 
are also evaluated at k = aH. The approximation scheme 
used to derive this result relies on assuming e is small as 
well as 77, and therefore the full strength of the exact 
power-law inflation result is not being used. 

We consider more general circumstances, by carrying 
out a general expansion about any power-law solution. 
Any general background model can be specified by giving 
the function H(<fi), or equivalently the value Hq — H(4>o) 
and the values of all the slow-roll parameters at horizon 
crossing for the mode of interest. We can then express 
H(cj)) in terms of the corresponding power-law expression 
(with the same values of H and e) and a series of cor- 
rection terms which depend on the higher-order slow-roll 
parameters: 
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The quantity we are interested in, V^[ 2 , is a functional 
of H((f)) and is therefore a function of the slow-roll param- 
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eters in the same form as they appear in the expression 
for H above: 

vU 2 =vU\e,a 2 ,a 3 ,...) (28) 
where we have defined the new set of parameters 
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These parameters give a measure of the 'distance' of the 
true model from the reference power-law model; they are 
all to be evaluated at k = aH so they are just numbers, 
not functions of <p. One should think of these as measur- 
ing the amount by which the derivatives of H at Hubble- 
radius- crossing fail to match those expected of the power- 
law solution given its amplitude and hrst derivative. For 
models close to power-law inflation these parameters will 

1/2 

be small and in such cases continuity demands that 
receives a small correction from the power-law result.^] 
This allows us to write the spectrum as an expansion in 
powers of small parameters; to first-order in (72,(73, .. . 
the expansion is 

V][ 2 (k) = [1 + A(e)a 2 + B(e)a 3 + ■ ■ -]V% 2 (k) (30) 

where P][ 2 {k) is the exact power-law solution Eq. (0) 
for the appropriate e and where everything on the right 
hand side is evaluated at k = aH. Here A(e), B(e), . . . 
are calculable functions giving, at each value of e, the 
coefficient in the expansion of each a n . In the limit of 
small e, A(e) must approach the value predicted by the 
Stewart-Lyth calculation 



A(e) 



C as e->0 



(31) 



where C = — 2 + ln2 + 7 ~ —0.73 is a numerical constant, 
7 being the Euler constant. 

In principle, to completely specify the background 
model we require values for an infinite number of slow-roll 
parameters, so even restricting the expansion Eq. j30| ) 
to first-order in each a n still leaves us with an infinite 
number of terms. However, as we have pointed out the 
perturbation spectrum depends on H(<fi) only for some 
interval of <fi around 4>q and so in a Taylor series expan- 
sion of H around 4>q we expect the higher-order terms 
(involving higher-order parameters <j n ) to be successively 
less significant. We therefore try a truncated form of the 
expansion Eq. (|30|), keeping only the first few terms. 

Unfortunately, we have not been able to determine 
A(e) analytically. Instead therefore we resort to a nu- 
merical computation. 



*For this to be strictly true, it is necessary that the func- 
tional P^ 2 depends only on a limited range of values of H((f>), 
which we have seen is the case here. 



III. NUMERICAL SOLUTION OF THE 
PERTURBATION EQUATIONS 

The perturbation equation is solved numerically. As an 
initial test, we check that the code can reproduce the an- 
alytical power-law inflation solution to high accuracy. In 
fact the asymptotic power spectrum output by the code 
differs from the analytical result by less than one part in 
10 5 . Fig. [I] shows a comparison between the output of 
our numerical code and the exact solution. 




* output of numerical code. 
— analytical solution 
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FIG. 1. A comparison of numerical and analytical solu- 
tions for perturbations generated by power-law inflation. 

It can be seen from Fig. || that before horizon crossing 
the perturbation amplitude decreases with time. At suf- 
ficiently early times the amplitude of the perturbation 
is large with respect to the background and the linear 
theory breaks down jl5| . We therefore do not expect the 
large initial values of the amplitude displayed in Fig. |] to 
be accurate. However, before horizon crossing the ampli- 
tude is fixed by the standard QFT normalization, Eq. (|8j), 
and is independent of the earlier behaviour, so since the 
perturbation becomes small before horizon crossing our 
results for the final amplitude using linear theory will be 
valid. 

Easther [[[9) has recently constructed another inflation- 
ary model for which the perturbation equations can be 
solved exactly; we have also successfully tested that our 
code can reproduce this solution. As it lacks many of the 
nice properties of power-law inflation, it does not appear 
promising for using as the basis of an expansion. 

Having thus confirmed the accuracy of the code, we can 
compare the numerical results with those obtained from 
analytical approximations based on the slow-roll approx- 
imation. As input to the code we use -ff (</>) as given by 
Eq. ( p7j ) and solve the perturbation equation for a range 
of values of e and 77 (or equivalently a 2) for the mode 
which crosses the horizon at <j) = cf>o. All higher-order er„ 
are set to zero. 

In Fig. |^ we display the relative differences between 
the numerical results and each of the four different ap- 
proximations we have considered here: the lowest-order 
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slow- roll expression, Eq. (g), the Stewart-Lyth calcula- 
tion, Eq. (p6[), the exact power-law solution for the cor- 
responding value of e, Eq. (]24|), and the expansion about 
the power-law solution given by Eq. (|30|). The region 
above the thick dashed curve in these contour plots has 
been excluded, because for these values of the parameters 
inflation ends soon after horizon crossing and the pertur- 
bation amplitude has not attained its asymptotic value. 
The shaded triangle in these plots represents the region in 
the e-rj plane favoured by current observations of large- 
scale structure and microwave background anisotropics 
pof , allowing for the uncertainty in cosmological param- 
eters. 

The final panel of Fig. || gives a comparison of the 
numerical results with the expansion Eq. ( |30|) and as such 
requires that we know the coefficient A(e). We have used 
the numerical code to evaluate this coefficient and the 
next-order coefficient B(e) and we find that for e < 0.3 
they are well fitted by the linear expressions 



A(e) ~ -0.730- 0.850e 
B{e) ~ 0.010 + 0.480e. 



(32) 
(33) 



As anticipated, we find the contribution from 03 to be 
significantly smaller than that from 02 ■ 

Before we compare the numerical results with the ap- 
proximations it is worth considering how representative 
our model is at each value of e and rj. In other words, 
how reasonable is it to neglect the effect of higher-order 
parameters? In order to get some idea of the impact of 
higher-order terms on the results, we have repeated the 
numerical calculations for the values of e and 77 shown in 
Fig. |^, using a different expansion for H as the input by 
instead expanding its logarithm: 
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but keeping only the first two terms. In this case we 
have specific non-zero values for all of the higher-order 
parameters in terms of e and er 2 j with 03 for example 
given by 

<Jz = 3a 2 e 1/2 . (35) 

We therefore have numerical results for two different 
slices through the parameter space. Comparing the re- 
sults for the two slices point-by-point in the e — rj plane 
we find that within the region shown in the contour plots 
of Fig. I the results differ by at most 2%. When Eq. (§|) 
is used to calculate the first-order correction arising from 
the contribution of er 3 , the relative difference between the 
two sets of results is found to be below 0.5% within this 
region. Of course it will always be possible to construct 
models where the higher-order terms are more significant 
than this, but in most cases the contributions arising 
from the parameters beyond 03 will be negligible. 



IV. GRAVITATIONAL WAVES 

We have also carried out calculations for the somewhat 
simpler case of gravitational waves, where the relevant 
equation, analogous to Eq. (||), is 



d 2 v\< 
dr 2 



k 



1 d 2 a 
a dr 2 



vk = . 



(36) 
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The corresponding power spectrum V g is expanded in 
the same manner as for scalar perturbations 

Vl' 2 {k) = [1 + A g (e)a 2 + B s (e)a 3 + ■ ■ -F g 1/2 (fc) , (37) 
~ 1/2 

where once again V g is the exact power-law inflation 
result and the coefficients are determined from numerical 
calculations to be 



B g (e) 



0.16e + 2.8e 2 
-l.le 2 . 



(38) 
(39) 



The absence of a constant term for A s (e) is in agree- 
ment with the Stewart-Lyth result for the gravitational 
wave spectrum. The differences between this approxi- 
mation and the numerical results are shown in Fig. |^, 
along with those for the other approximate expressions. 
Notice that in this case the power-law inflation result 
does better than the Stewart-Lyth result; this is because 
for gravitational waves, unlike density perturbations, the 
latter is simply the small parameter expansion of the for- 
mer. Note also that the results for gravitational waves 
are much more accurate than those for density perturba- 
tions. 



V. DISCUSSION 

In Fig. H the results based on the slow-roll approxima- 
tion and on the exact power-law expression are compared 
with the numerical results for the case of scalar perturba- 
tions. The lowest-order slow-roll result is seen to perform 
well only for very small values of the slow-roll parame- 
ters, with errors in excess of 10% within the observa- 
tionally favoured region. In contrast, the Stewart-Lyth 
result fares surprisingly well. The power-law expression, 
although exact along the line e = 77, rapidly decreases in 
accuracy away from this line. By construction, our new 
expansion about this exact solution shows the best agree- 
ment with the numerical results across the full range of 
parameters, though the improvement only sets in a con- 
siderable distance from the slow- roll limit. At small e 
and r\ the performance of our expansion is comparable 
with the Stewart-Lyth calculation. These results show 
that the assumption e<l, required for the Stewart- 
Lyth result to be applicable, doesn't break down until 
well outside the region constrained by observations. It 
seems then that for reasonable inflationary models the 
Stewart-Lyth calculation works extremely well. 
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Fig. H shows the comparisons between the numerical 
results and each of the approximations for the spectrum 
of gravitational waves. In this case the exact power-law 
expression is found to be extremely accurate. There is 
considerable improvement in expanding about this exact 
solution, but the high level of accuracy achieved is not 
likely to ever be required. 

As a check on the applicability of these results to gen- 
eral inflationary models, we have numerically evaluated 
the perturbation amplitude for a model with a steep poly- 
nomial potential, for the mode crossing the horizon 60 
e-foldings before the end of inflation. In this calculation 
the inputs to the numerical code are the potential and 
the initial conditions for the scalar field, the latter taken 
to be given by the inflationary attractor |f[s]| . We took 
the potential 



v{4>) = v 



(— 

\mp\ 



a 



10. 



(40) 



in order to be not too near the slow-roll limit; such poten- 
tials were discussed in Ref. [^2| . The slow-roll parameters 
take the values 



0.04 



r) = 0.03 . 



(41) 



The comparisons between the approximations and the 
numerical result show a relative error 



lowest order slow-roll 
Stewart-Lyth 
PLI 

new expansion 



5 x 10~ 3 , 

5 x 10" 4 , 

6 x 10" 3 , 
2 x 10~ 4 . 



(42) 



Note that for this case all approximation schemes fall 
within a 1% error. 

To recap, we have used the exact power-law result as 
the basis for an expansion which gives the perturbation 
spectrum very accurately for inflationary models which 
are 'close' to power-law models, in the sense that the pa- 
rameters <T ra are small. This expansion outperforms the 
Stewart-Lyth result for large e, but within the range of 
values for the slow-roll parameters allowed by cosmolog- 
ical observations we have shown that the Stewart-Lyth 
result is remarkably accurate. Our results indicate that 
the errors in both ours and the Stewart-Lyth calculation 
are due mostly to neglecting terms of order erf, with the 
contributions from higher-order parameters being rela- 
tively insignificant. We stress that these conclusions are 
only valid if the perturbation spectrum is sensitive only 
to the behaviour of H (</)) close to horizon crossing, but 
within the single scalar field paradigm we have been dis- 
cussing we expect this to be the norm and the result 
calculated for a specific potential is consistent with this. 
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FIG. 2. Contours indicating the relative differences between the numerical results and each of the approximate expressions 
for the amplitude of scalar perturbations. Dashed contours indicate an underestimate of the true result and solid ones an 
overestimate, by the percentage indicated. Above the thick dashed line the asymptotic regime is not reached before inflation 
ends, so no results are available. The shaded region indicates the parameters favoured by current observations. 
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FIG. 3. As Fig. ^ but for gravitational waves. 



9 



